Reservoir cross-over in entanglement dynamics 
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We study the effects of spontaneous emission on the entanglement dynamics of two qubits in- 
teracting with a common Lorentzian structured reservoir. We assume that the qubits are initially 
prepared in a Bell-like state. We focus on the strong coupling regime and study the entanglement 
dynamics for different regions of the spontaneous emission decay parameter. This investigation 
allows us to explore the cross-over between common and independent reservoirs in entanglement 
dynamics. 
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I. INTRODUCTION 

Since its discovery by Yu and Eberly in 2004, con- 
siderable interest has been devoted to the phenomenon 
of early-stage disentanglement, also called entanglement 
sudden death (ESD) [1]. Indeed not only is this phe- 
nomenon important from a fundamental point of view, 
but it is likely to play a role in many applications of quan- 
tum information theory and technology [2, For exam- 
ple, ESD might represent a threat to building a quantum 
computer. So far ESD has been predicted in many differ- 
ent theoretical systems such as pairs of qubits 0, la, 1,0, 
continuous variable system [8.], subset of multiple qubits 
and spin chains. It has also been observed experimentally 
in two different contexts • 

In general, a key factor determining the dynamics is 
the type of environment in which the system of interest is 
immersed. With just a pair of qubits, differences appear 
in the entanglement evolution, depending on the proper- 
ties of the environment. For example, when a bipartite 
system interacts with one reservoir (i.e. the two parts in- 
teract with the same reservoir), correlations between the 
two parts are created because of the reservoir-mediated 
interaction. These correlations cannot arise if each of the 
parts talks to only its own environment, and if there is 
no direct coupling present between them. On the other 
hand, if the correlation time of the bath is long, the envi- 
ronment keeps track of the dynamics of the small system, 
and revivals of entanglement may appear. These are the 
so-called memory effects, typical of non-Markovian dy- 
namics, appearing not only in common, but also in inde- 
pendent reservoirs. 

In Ref. [7] we have studied the exact dynamics of two 
qubits interacting with a common Lorentzian-structured, 
non-Markovian reservoir. Here we generalize these re- 
sults and take into account the spontaneous emission 
from the two qubits. This model describes, for exam- 
ple, the situation of two atoms strongly coupled to the 
same one-dimensional high-Q cavity, in which both the 
atoms can also independently emit a photon in any di- 
rection outside the cavity. The uni-dimensional high-Q 
cavity constitutes a common non-Markovian reservoir for 
the qubits, while the independent spontaneous emission 



of a photon in a flat continuum of modes can be seen as 
a consequence of the interaction with two independent 
Markovian reservoirs. 

The introduction of spontaneous emission for the two 
atoms allows us to study a rather realistic situation, char- 
acterized by different regimes of entanglement dynamics. 
In particular, we want to explore the cross-over between 
common and independent reservoir dynamics. 

In Ref. pjj] the generation of entanglement for two 
trapped ions coupled to a high-finesse cavity has been 
studied. In that case, not only the cavity losses, but 
also the entire atomic level structure (the ions are cou- 
pled to the cavity mode via a Raman scheme in a A- 
configuration) is taken into account. The aim of that 
work is to study realistic experimental conditions under 
which the collective Dicke model can be implemented in 
an ion-cavity QED context, so only states with one ex- 
citation in the atomic system are considered. Here we 
study entanglement dynamics when the atomic system is 
prepared in Bell-like states with one and two excitations. 

In the following sections we first introduce the mas- 
ter equation describing the system of interest (section 
mi) , then we study the entanglement dynamics for initial 
Bell-like states for different regions of parameters (sec- 
tion imp , and finally we interpret our results and make 
some conclusive remarks in section HVl 



II. THE MODEL 

We consider two two-level systems (qubits) interacting 
with the same leaky cavity in rotating-wave approxima- 
tion, and emitting independently outside the cavity due 
to spontaneous emission (see figure [1]). The qubits have 
the same transition frequency, and they are equally and 
resonantly coupled with the leaky cavity. In Ref. 0] 
we have studied the exact dynamics of two qubits in- 
teracting with the same Lorentzian-structured reservoir, 
i.e. a lossy resonator. In that work we found that, for a 
certain class of initial states, the entanglement dynam- 
ics exhibits regions of sudden death and resurrections. 
We have interpreted these results as a consequence of 
the memory effects of the non-Markovian environment 
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FIG. 1: (Color online) Diagrammatic representation of the 
model we consider. Two atoms couple to a cavity mode which 
leaks to the environment at a rate F. In addition, each atom 
can independently radiate directly to the environment at the 
rates 7a, 7b. 



C(t) 0. 




FIG. 2: (Color online) Concurrence as a function of scaled 
time and o? for two atoms prepared in the state and in 
the absence of spontaneous emission. (Taken from Ref. [2.].) 



and of the reservoir-mediated interaction between the 
qubits. For the sake of clarity, and as a reference point, 
we report in figure [2] the basic result for the concurrence. 
We solved the dynamics using the pseudomode approach 
The pseudomode master equation describes a 
coherent Tavis-Cummings type interaction between the 
qubits and the pseudomode, and the dissipative part in- 
volves the pseudomode only which leaks into a Marko- 
vian reservoir. For a high-Q cavity the pseudomode is 
basically identified as the cavity mode. 

Here we generalize such a master equation by adding 
phenomenologically the dissipative spontaneous emission 
part for the two qubits. Contrary to Ref. [3] we do 
not use the representation in terms of super-radiant 
and sub-radiant states, but we keep the natural basis 
{|00), |10), |01), |11)} for the two qubits, as the super- 
radiant basis has no advantage when the independent 
spontaneous emission processes break the symmetry of 
the system. Thus the master equation of the system in 



the interaction picture is 
d~p F 



dt 



i[H^ p] — —{a^ap + pa^a — 2apa^) 



-^{afa^p + pafa^-2a^paf), (1) 

where p is the density matrix of the qubits plus the cavity 
mode and 



(2) 



Here, and are, respectively, the Pauli raising and 
lowering operators for the atoms A and B, a and are 
the annihilation and creation operators for the cavity 
mode, Q is the cavity-qubits coupling constant, F is the 
cavity decay rate, and Ja/b are the spontaneous emis- 
sion rates for atoms A and B. For simplicity we consider 
7a = 7b =75 and the coupling Q to be the same for 
both qubits. 

We assume that the environment is at zero 
temperature, so there are at most two excita- 
tions in the total system. The basis we use 
to write the density matrix elements of p is 
{|000), |001), |002), |100), |101), |010), |011), |110)}, where 
the first and second digit indicate ground (0) and excited 
(1) state of qubit A and B respectively, and the third 
one indicates the number of excitations inside the cavity. 
We solve numerically the 64 differential equations for the 
density matrix elements, and then, to find the dynam- 
ics of the two-qubit system only, we trace out the cavity 
mode degree of freedom to find the reduced density ma- 
trix p. 

Once we have the reduced qubits density matrix we can 
derive the entanglement dynamics. To quan tify entan- 
glement we use the Wootters concurrence [ij , defined as 
C{t) = max{0, — a/As ~ a/AI}^ where {A^} are 

the eigenvalues of the matrix R = p{(jy®(jy)p*{(jy^(jy)^ 

with p* denoting the complex conjugate of p and cfy^^ 
are the Pauli matrices for atoms A and B. This quantity 
attains its maximum value of 1 for maximally entangled 
states and vanishes for separable states. 

We assume the qubits are prepared in Bell-like states. 
The evolution then drives such initial pure states into a 
mixed state having an "X" form 



/ a{t) w{t) \ 

h(t) z(t) 

z*(t) c{t) 

V ^*(t) d(t) J 



(3) 



For this class of states the concurrence assumes a simple 
analytic expression 



where 



C{t) = max{0,Ci(t),C2(t)}, 

C\{t) =2\w{t)\-2./m^, 
C2(t) =2\zit)\-2^ait)dit), 



(4) 

(5) 
(6) 



where coherences give a positive contribution to Ci{t) 
and C2{t) and hence to the concurrence, while the nega- 
tive parts involve populations only. Of course, any initial 
state in the "X" form of Eq. ([3]), pure, or mixed, would 
lead to the same expressions for the concurrence. For 
states not in this form we would have to perform a nu- 
merical evaluation of the concurrence. 



III. ENTANGLEMENT DYNAMICS 

We investigate entanglement dynamics for the follow- 
ing initial Bell-like states of the qubits: 

1$) =a|10) + e^^(l-a2)i/2|Q^^^ 

and 

1^) =a|00) + e^^(l-a2)^/2|ll). (8) 

With appropriate choices of and a these states include 
the usual four Bell states. 

Our aim is to understand the interplay between com- 
mon non-Markovian reservoir and independent Marko- 
vian reservoirs in entanglement dynamics. To do so we 
consider a strong coupling between the qubits and the 
cavity, and then we change the spontaneous emission de- 
cay rate 75. 

A Bell-like state with one excitation as in Eq. ([7j) never 
presents entanglement sudden death [7]. In the strong 
coupling regime and in absence of spontaneous emission 
the entanglement exhibits oscillations which are damped 
in time. The introduction of spontaneous emission from 
the two qubits leads to an additional damping in the 
oscillations. The result is that the larger the spontaneous 
emission rate is, the smaller the entanglement revivals 
are; for 75 > 3Q the oscillations are completely killed 
and the entanglement does not revive anymore. 

The dynamics of the Bell-like state involving two ex- 
citations, Eq. ([8|), is more interesting. First of all, let us 
recall the two limiting cases: no coupling with the cavity 
and no spontaneous emission. We obtain the Markovian 
independent reservoirs dynamics when setting the cou- 
pling between the qubits and the cavity to zero. In this 
case there are two different regions of the entanglement 
dynamics: for < 1/2 entanglement dies suddenly, and 
for > 1/2 entanglement decays exponentially. In the 
second limiting case, i.e. with cavity coupling and in ab- 
sence of spontaneous emission, the qubits interact exclu- 
sively with a common non-Markovian reservoir. In this 
case entanglement presents a much richer dynamics with 
oscillations for every value of and a series of dark peri- 
ods (ESD regions) and resurrections of entanglement for 
a'^ < 1/4, as shown in figure [2] [7]. 

Here we explore entanglement dynamics for different 
values of the spontaneous emission rate 75, fixing the 
strong coupling between the qubits and the cavity. In 
particular we set the cavity parameters to Q = O.270 
and r = VO.O570 where 70 is the decay parameter in 
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FIG. 3: (Color online) Concurrence as a function of scaled 
time and for two atoms prepared in state dS]). Here the 
spontaneous emission parameter 75 is ten times smaller than 
the cavity-qubits coupling constant Q. 



the case of an infinitely broad cavity. These parameters 
correspond to experimentally feasible conditions in the 
context of, for example, trapped ions [15] or circuit-QED 
\Tdl\ . The results do not depend on the relative phase 6. 

In figure [3] we present the entanglement dynamics for 
two qubits prepared in the state ([8]) with a weaker spon- 
taneous emission rate 75 = 1^/10. We notice that for 
short times the dynamics is not affected by spontaneous 
emission (compare with figure [2|). As time passes the 
effects of spontaneous emission become prominent: i.e. 
the oscillations get smaller in amplitude, and the ESD 
region dramatically increases until, for a time around 
t ^ 2/75 ^ IOO70, entanglement is dead irrespective of 
the value of a^. 

A clearer explanation of the dynamics can be found 
by looking at the evolution of the density matrix ele- 
ments. For a two-photon Bell-like state the concurrence 
is given by Ci(t), where the two-photon coherence gives 
a positive contribution, and the populations of the state 
with one excitation in one of the qubits give a negative 
contribution. Looking at the differential equations for 
the density matrix elements p, Eq. ([1]), we notice that 
the terms describing spontaneous emission cause a faster 
decay of two-photon coherence. On the other hand the 
population of the states |10) and |01) is less affected, since 
spontaneous emission not only adds a decay channel from 
|10)/|01) to |00) but also one from |11) to |10)/|01). As a 
consequence, at a certain time depending on the param- 
eter a^, the two-photon coherence will be smaller than 
the one-excitation populations and entanglement will be 
definitely lost. 

In the case of a strong spontaneous emission rate, 
75 = lOr^, the entanglement dynamics tends to be the 
same as the independent Markovian reservoirs case, as 
figure [4] shows. Even though the coupling with the cav- 
ity is strong, the entanglement does not revive and the 
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FIG. 4: (Color online) Concurrence as a function of scaled 
time and o? for two atoms prepared in state (|8]). Here the 
spontaneous emission parameter 75 is ten times larger than 
the cavity-qubits coupling constant Q. 




FIG. 5: (Color online) Concurrence as a function of scaled 
time and o? for two atoms prepared in state {|8]) for small 
(a) and long (b) time-scales. Here the spontaneous emission 
parameter 75 is equal to the cavity-qubits coupling constant 



spontaneous emission damps down any entanglement os- 
cillations completely. We also see that the entanglement 
vanishes exponentially for > 1/2 and dies suddenly 
for o? < 1/2. 

More delicate is the case of a comparable cavity-qubits 
coupling and spontaneous emission decay rate. Figure 
[5] shows the entanglement dynamics at different time- 
scales for 75 = Q. Compared to the case with no sponta- 
neous emission, the dynamics is quite different from the 
very beginning. Indeed the region of ESD appears to be 
strongly increased. However, later on the coherent inter- 
action with the cavity makes entanglement revive and, as 
time passes, revivals appear for a wider and wider range 
of the parameter, as shown in figure [5] (a). For longer 
times (figure [5] (b)), the range of a'^ over which we find 
multiple entanglement revivals shrinks, and then expands 
again and, in fact, we find that this pattern repeats many 
times. Clearly, as time passes the amplitude of the oscil- 
lations and revivals dramatically decreases, for example. 



for t ~ 150/70 the concurrence is of the order of 10~^. 
The details of the dynamics change a lot depending on 
the particular value of 75. However, for 75 of the or- 
der of 17, some general features of the dynamics can be 
identified. For example, for ^ and ^ 1 the 
dynamics of entanglement is mainly controlled by spon- 
taneous emission, resembling the independent Markovian 
reservoirs case. For intermediate values of a number of 
death and revival periods follow one after the other with 
decreasing intensity. This is a sign of the non-Markovian 
backaction of the cavity and of the cavity-mediated cou- 
pling between the qubits as well. In general the long- 
time picture looks like the superimposition of death and 
revivals on the independent Markovian reservoir dynam- 
ical pattern. 



IV. CONCLUSION 

We have studied the entanglement dynamics of two 
qubits interacting with a leaky cavity (shared non- 
Markovian reservoir) and emitting independently via 
spontaneous emission (independent Markovian reser- 
voirs). Starting with the exact master equation of Ref. 

describing the non-Markovian dynamics with the cav- 
ity, we have phenomenologically added two dissipative 
terms for the qubits. Depending on the ratio between 
the spontaneous emission parameter 75 and cavity-qubit 
coupling constant 17, different regimes in entanglement 
dynamics can be identified. When the spontaneous emis- 
sion rate is small, the entanglement dynamics is not very 
much affected for short times, but, as time passes, the re- 
gion of sudden death spreads out until for long times the 
entanglement is lost for every value of a'^. When 75 is 
much larger than the interaction with the cavity does 
not play any role. The two dissipative spontaneous emis- 
sion terms in Eq. ([1]) control the dynamics, preventing 
any possibility of revivals in entanglement. The inter- 
mediate region when 75 and Q are of the same order of 
magnitude shows a more complex behaviour exhibiting 
dynamical elements of both the two limiting cases. 

We emphasize that the method presented here can be 
applied to any general initial state of the qubits, and 
not only to Bell-like states. In particular, it is interest- 
ing to see the effects of mixedness of the initial state 
on the entanglement dynamics given that the dynam- 
ics depends dramatically on the state of preparation of 
the qubits. Here we just mention the case of an ex- 
tended Werner-like state with two excitations of the form 
p = r|^)(^| + (1 - r)I/4. This state is characterized by 
an entangled part, formed from the Bell-like state of Eq. 
([8]), and a maximally mixed part, identified with the unit 
matrix. Clearly, due to the mixedness of the state, the 
initial amount of entanglement is smaller than for the 
state 1^) alone. The presence of the mixed part causes 
a faster death of entanglement. However, for small val- 
ues of the 75 parameter the entanglement exhibits os- 
cillations and, as 75 increases, features of the reservoir 
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cross-over appear. 

In realistic experimental ion-cavity QED conditions it 
is not always possible to limit the losses of the system to 
cavity-losses, and often spontaneous emission needs to be 
taken into account. In our work we have demonstrated 
how the phenomenon of spontaneous emission comes into 
play in the entanglement dynamics of strongly interacting 
cavity QED systems. 
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